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Abstract 


An Engel series is a sum of the reciprocals of an increasing sequence of 
positive integers, which is such that each term is divisible by the previous 
one. Here we consider a particular class of Engel series, for which each 
term of the sequence is divisible by the square of the preceding one, and 
find an explicit expression for the continued fraction expansion of the sum 
of a generic series of this kind. As a special case, this inclndes certain series 
whose continued fraction expansion was found by Shallit. A family of 
examples generated by nonlinear recurrences with the Laurent property is 
considered in detail, along with some associated transcendental numbers. 

Keywords: continued fraction, nonlinear recurrence, transcendental num¬ 
ber, Laurent property. 
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1 Introduction 

Given a sequence of positive integers (x„), which is such that Xn\xn+i for all n, 
the sum of the reciprocals is the Engel series 



( 1 . 1 ) 


where yi = xi and t/n+i = Xn+ijxn for n > 1. (It should be assumed that 
{Xn) is eventually increasing, which guarantees the convergence of the sum. A 
brief introduction to Engel series can be found in [^.) In recent work |12j . 
we considered some particular series of this kind that are generated by certain 
nonlinear recurrences of second order, and are such that the sequences {xn) 
and (j/ri) appear interlaced in the continued fraction expansion. Here we start 
from a sequence with the stronger property that x'^\xn+i- Any initial Is can be 
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ignored, but in what follows it will be convenient to start with xi = 1 and take 
X 2 > 2, which implies Xn > 2^ for n > 2, and we may write 

n 

Xi = l, zf " (1.2) 

i=2 


for some sequence of positive integers (z„) with Z 2 > 2. 
Engel series is 




, 23-2 23-3 




3=2 


The corresponding 


(1.3) 


Our main result will be to show that in the generic case, when Z 2 > 3 and > 2 
for n > 3, the continued fraction expansion of S has a universal structure, which 
we present explicitly. 

There is a precedent for these results in the work of Shallit, who first found 
the continued fraction expansion of the sum 


E 





(1.4) 


in m. for integer rt > 3 (with u = 2 being a degenerate case), and went on [18] 
to obtain the continued fraction for the more general series 


OO 


E 


1 


(1.5) 


where (cfc) is a sequence of positive integers with some non-negative N such 
that dn ■= Cn+i — 2c„ > 0 for all n > N. If we set Z 2 = u'^° and Zj = it'^3-3 for 
j > 3 in (II.3L and assume d„ > 0 for all n > 0, then ^ — 1 coincides with ( 11 . 51 ) . 


1.1 Outline of the paper 

In the next section we prove the main result, namely the expression for the 
continued fraction expansion of a generic series of the form (11.31) . Section 3 is 
devoted to an infinite family of examples of series of this type, which are gen¬ 
erated by nonlinear recurrences with the Laurent property. For such nonlinear 
recurrence sequences, we show that the sum of the corresponding series (IE3 
is a transcendental number. In the fourth section we consider the continued 
fractions obtained from certain degenerate cases, when either Z 2 = 2 or = I 
for n > 3, and some particular examples of these degenerate cases are examined 
in more detail. The final section contains some conclusions. 
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2 Continued fractions 

We use the notation 


[tto; 01 , 02 ,03, . . . ,Oj, ...] = Oo + 


Ol 


02 


03 


aj + ■ 


for continued fractions, and for the nth convergent of the continued fraction 
[oo; Ol, 02 ,...] we have 

Pn , 1 

— — Oo; Ol, 02 ,..., o„j, 

Qn 


where the numerators Pn and denominators qn are given in terms of the coeffi¬ 
cients ttj according to the matrix identity 


/oq 1\/oi 1\ / Qji 1 \ _ / Pn Pn — 1 A /"O 1 'l 

Vl Qn-l ) ■ ^ ’ 

For what follows, it will also be convenient to note the identity obtained by 
taking the determinant of each side of (EU, that is 

PnQn-l - Pn-iq-n = ■ (2.2) 


For other basic results on continued fractions, the reader is referred to [3]. 

To begin with, it is convenient to treat the factors Z 2 , zs, ■ ■ ■ of the sequence 
(Xn) as variables. For the first few partial sums we find the continued fraction 
expansions 


5i = 1 , ^2 


Z2 + 1 
Z 2 


[i;^ 2 ], S 3 


Z2Z3 + Z2Z3 + 1 
Z|Z 3 


[1;Z2-1,1, 2:3-1, ^2], 


where the nth partial sum of is denoted Sn- In general it is straightforward 
to write Sn as a fraction, that is 


5. = 


= l 11A;=: 


J 

k ^2 


‘e_ 2 J- 






• • Zn 


(2.3) 


where the denominator is Xn as given in (13); but the continued fraction ex¬ 
pansion of the nth partial sum is best described recursively. 

The basic pattern can be seen by looking at the continued fraction for the 
fourth partial sum, which is 


5'4 = [1; 2:2 - 1,1, 2:3 - 1, 2:2, 2:4 - 1,1, 02 - 1, 2:3 - 1,1, 2:2 - 1]. 

Observe that the first five coefficients are the same as those of S 3 , followed by 
Z 4 — 1 , 1 , and then four more coefficients which almost coincide with the last 
four in S 3 in reverse order, except that there is 2:2 — 1 in place of the final Z 2 in 
S 3 . This pattern persists, as described by the following 
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Proposition 


2.1. Given the initial set of coefficients 

r (3) (3) (3) (3) (3)i -I 

K ;«! ,02 :«3 J = [1;^2 - 1,1, 2:3 


1 ,^ 2 ] 


of the continued fraction expansion of S 3 , of length £3 = 4, define a sequence 
of sets of coefficients with £„= 3 ■ 2"“^ — 1 for n = 3,4,5,... 

recursively according to 




for j = 0, - 1, 

— r — 1 — 1 — 1 

— Zn +1 r, 0^„ + l — -L, 0£„+2 ~ 


and 


for j = in + 3,..., 2£„. 


r/ien t/ie nth partial sum of the series has the continued fraction expansion 


'S'n = [ao"\aJ"\...,4"Li]- (2.4) 

Proof: This can be done similarly to the proof in m, but we prefer to use 
matrix computations, in the same vein as m- The case n = 3 is easily verified 
directly. Proceeding by induction, suppose that the continued fraction expan¬ 
sion of Sn is given by (1^ . with coefficients for _) = 0 ,... in — I with 
in = 3 ■ 2”“^ — 1 defined according the prescription above, and denote the nu¬ 
merators and denominators of the convergents by pj and qj respectively; so the 
final convergent gives 

c _ _ 

— 5 Q^ri — ^ — ■ 

Qin-l 


Then for the next finite continued fraction defined by this recursive procedure 
there are a total of £n+i = ^in 4-1=3- 2”“^ — 1 coefficients, that is 
for j = 0,... in+i — 1, and for the convergents we use Pj,qj to denote numera¬ 
tors/denominators, respectively. So by (j2.1[) we have 


W„+1-1 

W„+l-2 

9^n + l-l 

9C+1-2 


= M„ 


^n+l 1 1 

1 0 


1 1 
1 0 




where 




1 0 


af"Ll 1 ^ ^ Wn-l Wn-2 \ 

1 0 / ~ V J 


(2.5) 


and 



_ 1 


„(”) _ 1 
1 


1 

0 

1 

0 




9^n-2 



Wn-l ~ ®n-l + ®n-2 ” W„-2 

w„-2 - gc-2 


-1 
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with T denoting transpose. Thus the equation (12.51) simplifies to yield 


f W„+1-1 Pl„+i-2 \^( + 1 2:n+lW„-lA„ - 1 \ 

V 9^„+l-l %n+l-2 / ~ V ^n+lQi^-i Z„+i®„_iA„ - 1 / ’ 

with A„ = pt^-i — %„-i, where we have used the fact that 


det M„ 


W„-i Wn-2 
®„-l 9^„-2 


= -1 


by (12.21) . since is odd. Hence we have 


W„+i-l — 2:n+l9f„-lPC-i + 1> 


-1 — — X — ^n+l5 


SO that ^ 

‘S'n+l — Sn H 

^n+1 

which is the required result. 


Pf.n -1 


1 




P^n + l-l 

%n+l-l 


□ 


Remark 2.2. Note that, mutatis mutandis, both the recursive structure of the 
partial sums and the above inductive proof hold for the partial sums of an Engel 
series cny if for some positive integer no, the sequence (xn) satisfies the weaker 
condition that x‘^\xn+i for n> no only. This is the analogue of the fact that for 
the series ms) in m, Cn+i — 2cn > 0 need only hold for n > N, for some N. 

The finite continued fraction expansions of the partial sums immediately yield 
the continued fraction for the full series 113, at least for a generic choice of 
factors Z 2 , 2 : 3 ,... of the sequence (a;„). 

Theorem 2.3. For integer factors Z 2 > S and Zn > 2 for all n>3, the Engel 
series {H has the continued fraction expansion 

S = [oo;01,02,... ,Oj,...] = [l;z2 - 1 ,1,2:3 - 1,2:2,2:4 - 1 , 1 ,...] ( 2 . 6 ) 

where the coefficients are given by aj = lim„_>.oo o^"^. 

Proof: The result follows from taking the limit n —>■ 00 in (1^ . provided that 
none of the coefhcients in the finite continued fractions vanish. To see that 
this is so, note that only 1 ,Z 2 ~ 1 , 2:2 and Z 3 — 1 appear as coefhcients in the 
continued fraction for S 3 , and all of these are non-zero with the above conditions 
on the Zj. At each step of the recursion in Proposition 12.II only z„+i — 1 > 1 
and a !'^'^2 — ~ 1 potentially new coefhcients, so we must check that 

— 1 cannot vanish. For n = 3 the last coefficient is a!f^ = Z 2 , so in 54 
this gives a!^'^ = Z 2 — 1, while for n > 4 we have = ^2 — 1, so 

af"+ 2 ^ = ^2 — 2, which is why we require Z 2 > 3. Thus the only numbers that 
appear as coefhcients in the continued fraction expansion of S are 1 , 22 , 2:2 — 2 
and Zj — 1 for j > 2, and none of these are zero. □ 
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3 Nonlinear recurrence sequences 

Among nonlinear recurrences of the form 

^n+N — f (^n+1; ■ • ■ ; ^n+N—l ); (^* ^) 

where / is a polynomial in A^ — 1 variables, there is a multitude of examples 
which surprisingly generate integer sequences. In a wide variety of cases, the 
recurrence m has the Laurent property: for certain special choices of /, all 
of the iterates belong to the ring ..., as a consequence, if all 

the initial values are 1 (or ± 1 ), then each term of the sequence is an integer. 
Such sequences were popularized by Gale mis], and subsequently Fomin and 
Zelevinsky found a useful technique - the Caterpillar Lemma [7] - which can be 
used to prove the Laurent property in many cases, i.e. for recurrences coming 
from cluster algebras [ 6 ] or in the more general setting of Laurent Phenomenon 
(LP) algebras [H] . 

In [TUI E] we classified recurrences of second order, of the form 

Xn+2 Xn = /(x„+i). (3.2) 

For the Laurent property to hold, the recurrence (13.21) must belong to one of 
three classes, depending on the form of /: (i) /(O) ^ 0, in which case one can 
apply the framework of cluster algebras (when / is a binomial) or LP algebras 
(when it is not); (ii) /(O) = 0, /'(O) ^ 0; (iii) /(O) = /'(O) = 0. In the first two 
classes there are additional requirements on /, but in the third class one can 
take /(x) = x^F(x) with arbitrary F G Z[x]. 

In [TU] we considered the case that /(x) = x^F(x), where F has positive in¬ 
teger coefficients with F'(O) = 1, and obtained the continued fraction expansion 
of the sum order to obtain an Engel series of the form (II.3|) . we 

should instead choose F so that (ESI) becomes 

Xn+2Xn = X,^_^-^G(^Xn-\-l\ 


where 


di>3, G(x)gZ>o[x], degG = d 2 > 0 , G(0) 7 ^ 0, G(l) > 3. (3.4) 


From Xn+ 2 /xf^^i = ■x‘^l^^G{xn+i), we see by induction that, starting with 

the initial values xq = Xi = 1 , (x„) is a sequence of positive integers such that 
x^|x„+i with X 2 = G(l) > 3; hence also = x„/x^_]^ > 3 for n > 2. 

Example 3.1. Taking di = 3, and G(x) = 3 for all x (so d 2 = 0), the re¬ 
currence fSl) becomes Xn+ 2 Xn = 3x^_|_]^, which generates the sequence begin¬ 
ning 1,1, 3,81, 531441, 5559060566555523,.... In this case the recurrence can be 
solved explicitly to yield 


5-75/3-^75^ 5-b75 /3-75\ 

“ 10 2 j 10 2 j 


n 
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The sum of the reciprocals is 


S'—1 + - + ^ + H-— [1; 2,1,8,3,80,1,2,8,1,2,19682,...]. 

Remark 3.2. In the above example, S — 1 is a sum of the type considered 
in m, with u = 3 and Cn = Sn- 2 - More generally, any recurrence of the form 
Xn+2Xn = di > 3 and xq = xi = 1 generates a sum of this type. 

For most choices of G it is not possible to give the general solution of the non¬ 
linear recurrence (13.311 in closed form. Nevertheless, one can adapt the methods 
of Aho and Sloane [5] to write a formula giving precise asymptotic information. 
By rewriting (13.31) in terms of logarithms we find that A„ = log satisfies 

A„+i - (di d 2 )A„ A„_i = log c -b with a„ = log ( -^ j , (3.5) 

where G[x) = cx'^'^ -I- lower order. Since = log(l + 0{xf^)) = 0{xf^) as 
n —>■ oo, the leading order behaviour of A„ is determined by the linear expression 
on the left-hand side of (13.5L which has the characteristic equation — {di + 
d 2 )X -1-1 = 0, with largest root 

X=^1±±±^/Mi±WEI^2. (3.6) 

The next two statements are equivalent to analogous formulae obtained for the 
sequences considered in [12]. 

Proposition 3.3. For the initial conditions Xq = Xi = 1, the logarithm A„ = 
log Xn of each term of the sequence satisfying i3.3\) is given by the formula 


A„, = 


(1 - A-i)A’^ - (1 - A)A- 


A-A 


-1 


n —1 


— 1 ) log C 




yn—k _ 


A-A 


-1 


(3.7) 


where otk is defined as in and A as in Eg). 

Corollary 3.4. To leading order, the asymptotic approximation of the logarithm 
A„ is given by 


A 


n 


CA" 


(3.8) 


where 


G = 


di d2 — 2yA — A ^ 


1 - A 


-1 


logc -b 


.. OO 


A-A 


k=l 


and for the terms of the sequence 


c d'i+<i 2-‘2 exp(CA"). 


The asymptotic behaviour of these nonlinear recurrence sequences is enough 
to show that the sum of the corresponding Engel series is transcendental. 
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Theorem 3.5. Suppose that the sequence (x„) with initial values xq = xi = 1 
is generated by the recurrence iS.Sl) for some G satisfying the conditions jS.fl )- 
Then the sum S in is a transcendental number. 


Proof: This is essentially identical to the proof of Theorem 4 in m, so here 
we only sketch the argument. Recall that Roth’s theorem says that if a is an 
irrational algebraic number then for an arbitrary fixed i5 > 0 there are only 
finitely many rational approximations p/q for which 


P 

a - 


q 


1 

q2+S • 


(3.9) 


The number S has an infinite continued fraction expansion, so it is irrational. 
From the asymptotics (13.811 it follows that for any e > 0 the growth condition 


Xn+l ^ X.j^ 


(3.10) 


holds for all sufficiently large n. By making a comparison with a geometric sum, 
this gives 

S_PJ^ 

qtn-l 


= E 


< 


—€—€ 




for any e' > 0 and n large enough. So if e and e' are chosen such that X — e — e' = 
2 + 5 > 2, then a = S has infinitely many rational approximations satisfying 
(j3.9L and hence must be transcendental. □ 


Example 3.6. Taking di = 3, and G{x) = 2x + 1 for all x (so d 2 = 1), the 
recurrence fSJ) becomes 

^n+ 2 ^n ~ T 1)5 (^-H) 


which generates the sequence 

1,1,3,189,852910317,5599917937724687764238078261637795,... 
with leading order asymptotics 

gC(2+V3)" 

Csi 0.107812043. 

The sum of the reciprocals is the transcendental number 

S = [1; 2,1,20,3,23876,1,2,20,1,2,7697947188058154,...] « 1.3386243. 

Engel series of the form can also be generated by nonlinear recurrences 
of higher order. For instance, one can take a recurrence of third order. 


V2 


X„+ 3 X„ = X^Yl^n+ 2 ^(^"+l,^n+ 2 ), with Cl > 1 , 63 > 2 


(3.12) 











where the polynomial H{X,Y) G Z>o[X, F] is not divisible by either of its 
arguments. It is straightforward to show that the Laurent property holds for 
this recurrence, and from X„+ 3 /X^_j .2 = -^n+ 2 ), h 

is easy to see by induction that the initial values Xq = Xi = X 2 = 1 generate 
an integer sequence with X^\Xn+i for all n > 0. Thus the sum of reciprocals 
starting from the index 2 , that is 




t=2 




(3.13) 


is an Engel series of the form (11.31) . Note that the condition > 3 should 

be imposed, in order for Theorem 12.31 to apply to this series. 

A particular class of recurrences of the form (13.121) can be obtained by fac¬ 
torizing the terms of a sequence satisfying (1531) as Xn = XnXn+ 1 , which lifts 
the second order recurrence to 


A„+3A„ = (A„+iA„+2)^i-1G(X„+iA„+2). (3.14) 

For a generic polynomial iL(X„+i, A„_|_ 2 ) on the right-hand-side of (13.121) . it is 
not immediately obvious which term will be dominant as n —>■ 00 , but in the 
special case (13.141) the same techniques as for the second order recurrence can 
be applied directly, to show that the leading order asymptotics is log A„ ^ G'A" 
for some G' > 0, where A is given by (13.6|) . This means that Theorem [33] applies 
to the series (I3.13p as well. 

Example 3.7. Setting Xn = XnXn+i in iS.lll) gives the recurrence 

A„+3A„ = A2^iA2^2(2A„+iA„+2 + 1), (3.15) 

which generates the sequence beginning 


1,1,1,3,63,13538259,413636490314204194515563505,.... 

To leading order, log ~ G'(2 -|- v^)", with C ~ 0.0227833. The sum of the 
reciprocals in h3.13\) is the transcendental number 

S' = [1; 2,1,6,3,3410,1,2,6,1,2,2256800700104,...] « 1.3492064. 

For other examples of transcendental numbers whose complete continued 
fraction expansion is known, see [4] and references. 


4 Degenerate cases 

If either Z 2 = 2 or = 1 for some n > 3, then one of the coefficients in the 
continued fraction becomes zero, and Theorem l2.3l is no longer valid. To obtain 
a continued fraction with non-zero coefficients, one can use the replacement rule 
[..., a, 0, 6 ,...] —>■ [..., a -I- &,...] (see Proposition 3 in [T5]) to remove the zero. 
Each such replacement, decreases the length of a finite continued fraction by 
two, so in degenerate cases the length of the continued fraction expansion of 
is typically shorter than the generic value £n = ^ ■ 2"“^ — 1 . Here we present 
the expansion for two particular degenerate cases, omitting details of the proof. 
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4.1 The case Z 2 = 2 


For generic values of the factors Zj , the sequence of lengths in of partial sums 
begins 1, 2, 5,11, 23 for n = 1, 2, 3,4, 5. When Z 2 = 2 and Zj > 2 for all j > 3, 
the first few continued fractions for the partial sums of S are S'! = 1, 

S 2 = [1-,2], ,S3 = [1;1, 1,^3-1,2], 54 = [1;1,1,^3-1,2,Z4-1,1,1,^3-1,2], 

which are of the same length as in the generic case, except for S4 being of 
length 10, since at the end , 1,1] —>■ , 2]. The first zero appears in S' 5 , 

which contains a single coefficient Z 2 — 2 , so removing this and making the final 
replacement , 1 , 1 ] —, 2 ], results in the length being 20 : 

S5 = [1;1, 1,2:3-1,2,Z4-1, 1,1,Z3-1,1, 1, 2:5-1, 2, 2:3-1,1,1, 24-1, 2, Z3-1, 2]. 

Thereafter the pattern continues with the continued fraction doubling in length 
at each step, as described by the following 

Theorem 4.1. When Z 2 = 2 and zj 7 ^ 1 for j > 3, the Engel series hl.3\) has 
the continued fraction expansion 


S = [oo; 01,02, ...,aj,...] = [ 1 ; 1 ,1,2:3 - 1 , 2 ,2:4 - 1 , 1 ,...] 

with coefficients given by Uj = lim„^oo 0^"^ where sequence of 

sets of coefficients of the finite continued fractions for partial sums Sn, of length 
in = 5 ■ 2 "“^ for n = 4, 5, .. is defined by starting from 

[o[)'‘^ 0^"^^ ..., 0^"^^] = [1; 1,1, 2:3 - 1, 2, Z4 - 1,1,1, 2:3 - 1,2], 


and obtaining subsequent coefficients according to 


(ra+l) (n) 

0 } >= a) > 


for j = 0 ,...,in- 2 , 


(n+l) _ 

I — -L, 


ol"+') = 1, 




(n+l) _ (n) 

— “ 2^„-4 + 1 


= 2 . 


n + l ~ Zn+1 — 1, 

for j = in + 2 ,, 2 in — 2 , and 

In order to obtain a sequence (Xn) of this degenerate type from a second 
order recurrence of the form (liOll . the conditions (ITl) should be modified so 
that G(l) = 2, which requires that G(x) = x‘^^ + 1 for some non-negative integer 
d 2 . So the recurrence becomes 


Xn+2 Xn = X^+iix^+l + 1), (4-1) 

with di > 3, d 2 > 0. The results of Proposition 13.31 Corollary [TH and Theorem 
13.51 all apply without alteration to sequences obtained from 6ID. 

Example 4.2. Taking di = 3, and G{x) = a; -|- 1 for all x, the recurrence I3.3\) 
becomes Xn+ 2 Xn = Xn+i{xn+i + 1), which generates the sequence 


1,1,2,24,172800,37150633525248000000,... 

with asymptotics Xn ^ gC{2+V3) ^ (j ~ 0.06224548. The sum of the reciprocals 
is the transcendental number 


S = [1; 1,1,5,2,299,1,1,5,1,1,1244167199,2,5,1,1,299,...] r; 1.54167245. 
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4.2 The case Zj = 1 for j > 3 

If we set Z 2 = u and all other factors Zj = 1 then for n > 2; the 

expansion isrm is no longer valid because each coefficient Zj — 1 becomes zero 
for j > 3. In that case, the sum of the reciprocals is S' = 1 + 
that S — 1 coincides with m- The continued fractions for the partial sums 
were first obtained in |16j , and a nonrecursive description was given in |17| . The 
sequence of lengths begins 1,2,3,5,9,17, with = 2"“^ + 1 for n > 3, and the 
full continued fraction is 

S = [l;u — l,u + 2,u,u,u — 2,u,u + 2,u,u — 2,u + 2,u,u — 2,u,u,u + 2,u,...] 

for u ^ 2. The only numbers that appear as coefficients in this continued 
fraction are 1 , u — 2 , u — 1 , m, tt + 2 . 

However, the case u = 2 is special, since some of the coefficients in the above 
expansion become zero. The sequence of lengths of partial sums starts with 
1,2,3,5,7,11, and = 2"“^ + 3 for n > 4. Only the numbers 1,2,4 ,6 appear as 
coefficients in the continued fraction for the series, which is 

S = [1; 1,4, 2,4,4, 6 ,4, 2,4, 6 , 2,4, 6 ,4,4, 2,4, 6 ,...]. 

The argument used to prove Theorem 13.51 based on Roth’s theorem, does 
not apply to the partial sums of the series m, since the sequence (x„) does not 
grow fast enough: in contrast to (I3.10L it satisfies the recurrence Xn+i = 
However, a direct proof of transcendence of (d, valid for all integers u > 2 , 
was first given by Kempner [13]; various alternative proofs are collected in [Tj. 


5 Conclusions 

We have found the continued fraction expansion for an Engel series of the special 
type d- In some cases, coming from nonlinear recurrence sequences, it has 
been shown that this produces transcendental numbers. We expect that the 
sum d should be transcendental for any choice of the factors Z 2 > 2, and 
Zj > 1 for j > 3. However, we do not know of a simple way to prove this in 
general. 

Acknowledgments: This work is supported by Fellowship EP/M004333/1 
from the EPSRC. The original inspiration came from Paul Hanna’s observations 
concerning the nonlinear recurrence sequences described in |12j . which were 
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